How to Directly Measure Kondo Cloud's Length 
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We propose a method to directly measure, by electrical means, the Kondo screening cloud formed 
by an Anderson impurity coupled to semi-infinite quantum wires, on which an electrostatic gate 
voltage is applied at distance L from the impurity. We show that the Kondo cloud, and hence the 
Kondo temperature and the electron conductance through the impurity, are affected by the gate 
voltage, as L decreases below the Kondo cloud length. Based on this behavior, the cloud length 
can be experimentally identified by changing L with a keyboard type of gate voltages or tuning the 
coupling strength between the impurity and the wires. 

PACS numbers: 72.15.Qm,73.63.Kv,72.10.Fk 



Introduction. — The Kondo effect is a central many- 
body problem of condensed matter physics [H, • It in- 
volves a spin singlet, formed by the spin-spin interaction 
between a magnetic impurity and surrounding conduc- 
tion electrons. Deeper understanding of the effect has 
been achieved by using a quantum dot, that hosts a mag- 
netic impurity spin, under systematic control 

Although the Kondo effect is well known and widely 
studied, its spatial features still remain to be addressed. 
The Kondo spin singlet is formed below a certain energy 
scale, the Kondo temperature Tk- This implies that the 
singlet is formed over a spatially extended conduction- 
electron region of length scale £k = hvp/ (fcgTR-), where 
vf is the Fermi velocity. The region is called the Kondo 
screening cloud. When Tk ~ IK and vf ~ 10 5 -10 6 m/s, 
the cloud has the size of micrometers. There have been 
several proposals for ways to detect the cloud. 

Despite of the proposals, there has been no conclu- 
sive measurement supporting the existence of the Kondo 
cloud 0, [2l[ . The difficulty to detect the cloud arises 
because it is essentially a spin cloud. More than that, the 
conduction spin forming the Kondo spin singlet exhibits 
quantum fluctuations, resulting in zero averaged spin. A 
way of detecting the cloud is to study the spin-spin cor- 
relation (71-Hoj between the impurity and the conduction 
electrons. However it requires measurements of spin dy- 
namics of time scale S/^T/f). STM studies probing 
local density of states may be also useful [Tl - 13 1 . Recent 
STM measurements |22h24| have progress, showing the 
signatures of the Kondo effect in the region away from 
a magnetic impurity, whose spatial extension is however 
much shorter than Another direction is to study a 
magnetic im pur ity in a finite-size system, using persis- 
tent current 14 1 or electron conductance 15 - IH . In this 
case, because the cloud cannot extend beyond the finite 
size, the Kondo effect is strongly affected, and suppressed 
when the size is shorter than £k- There has been no ex- 
perimental report in this direction. 

In this work, we propose a new way of detecting the 
Kondo cloud, based on the intuition that a change of con- 



duction electrons inside the cloud will affect the Kondo 
effect. We consider a Kondo impurity formed in a quan- 
tum dot (at x = 0) coupled to two semi-infinite ballistic 
quantum wires with electron tunneling amplitude twD 
(see Fig. [TJ). Electrostatic gate voltages V g are applied 
to the wires at distance L from the dot (in \x\ > L), 
and modify indirectly the local density of states p(e) of 
conduction electrons nearby the dot (see Fig. [2J ; our ap- 
proach is applicable also when V g is applied to only one 
wire. Wc find that when L 3> V g does not affect 
the cloud. However, when L -C £k, the cloud, hence the 
Kondo temperature Tk and the electron conductance G 
of the dot, are sensitive to V g . The crossover between 
the two regimes occurs at L ks £k- By measuring G or 
Tk with varying L or ty^D (see Figs. [3] and [4]), one can 
detect the crossover and hence £k- We remark that the 
cloud is extensible to \x\ > L, hence, not suppressed even 
for L <C £k in our approach; this is in contrast to the 
case 13- 17 1 where the cloud is limited by the finite sys- 
tem size. We analyze the system using the poor man 
scalin g |25f , t he numerical renormalization grou p (NRG) 
study [26h28| . and the Fermi liquid theory [29(. 

The setup. — We describe the wires by the tight- 




FIG. 1: (Color Online) Setup for detecting the Kondo cloud. 
A quantum dot located at x = hosts a magnetic impurity 
spin. It couples to quantum wires along x axis, with electron 
tunneling amplitude twD- Gate voltage V g is applied at dis- 
tance L from the dot (in \x\ > L). The Kondo effect becomes 
sensitive to V g , as L decreases below the cloud length. 
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binding Hamiltonian with sites j's in wire i = l,r, 

oo 

H W = E foray* + M C Iia C *(j + l)<x +H.C.)] 

i=J,r j=l cr=f,4. 
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where c| JCT creates an electron with spin a and energy 



eo at site j in wire i, riijo = c\- a Cij a , arid i is the hop- 
ping energy. The last term describes the gate voltage 
V g applied to the two wires in \x\ > L = Na, where a 
is the lattice spacing. V g (x) changes at x = L abruptly 
over the length shorter than the Fermi wave length. The 
dot Hamiltonian Hr> is modeled by the Anderson im- 
purity doj, H D = Y^(r=t,i e ddkd a + Un^ridi, where d\ 
creates an electron with energy e d and spin a in the dot, 
n da = do-dp, and U is the electron repulsive interaction of 

the dot. H T = -t W D J2 t Z)<r=t,4.( c lio- d <>-+ H - c -) describes 
electron tunneling between the wire and the dot. 

The dot is occupied by a single electron in the Coulomb 
blockade regime of e d < e F and T(e F ) <C — e d + ef, U + 
£d — £F) where T(e) = 2ir\t\y d\ 2 p(e) is the hybridization 
function between the dot and the wires (i.e., the level 
broadening) , p(e) is the local density of states at energy 
e in the sites j = 1 (the neighboring sites of the dot) of 
the wires, and ep is the Fermi energy. In this regime, 
the total Hamiltonian of the setup is transformed by the 
Schrieffer-Wolf transformation into @, El 

H = H D + H w + H T 2± Js ■ S + n u<7 + H w- (2) 



Here, the Kondo impurity spin, S = a< ^a^ao'd-a 1 /2, 
couples with the spin of the neighboring conduction elec- 
trons, S = Y.o.a'( C \la + c llo)Vcrcr' (Cllcr> + C rla <)/2 7 with 

strength J = 2t 2 wr ,[-l/ {e d - e F ) + 1/(17 + e d - e F )}. 
The second term describes the potential scattering with 



strength V = t 2 WD [-l/{e d - e F ) - l/(U + e d - e F )}/2. 

Local density of states. — We show how the gate volt age 
V g changes the local density of states p(e) of conduction 
electrons near the impurity (at sites j = 1). We cal- 
culate p(e), by matching the wavcfunctions of the wire 
Hamiltonian Hw between j = N and N + 1, 



Pit) 



sin(g'o) sin 2 (fca) 



nat[sm 2 (ka) + ^f- sin[fc(iV + l)a] sin(fcAa)] 



(3) 

k and q are the wavevectors in |x| < L and \x\ > 
L, respectively, satisfying e = eo — 2icos(fca) = eo — 
eV g — 2tcos(qa). For finite L, the hybridization function 
T(e) = 2irt WD p(e) shows resonances with level spacing 
A = irhv F /L; see Fig. [2] The oscillation amplitude of 
T(e) is proportional to V g . The average value of T(e) 
around e F is denoted as Too. 




FIG. 2: (Color Online) Hybridization function T(e) — 
2jvtwDP( € ) f° r infinite (blue dashed curve) and finite L (green 
solid); we choose e /2i = 0.931 and eV g /2t = 0.125. For fi- 
nite L, F(e) shows resonances with spacing A = nhvF/L; 6f is 
chosen to be located at a resonance center (see the red dashed 
line) . When L is so large that V g is applied outside the Kondo 
cloud, the Kondo temperature Tk = Tkoo is determined by 
Too and independent of V g . On the other hand, Tk depends 
on V g and L, when L is smaller than the cloud size £k- When 
L is much smaller than £k, Tk is determined by To = F(eF)- 



We sketch how the change of p(e) by V g affects the 
Kondo effect in different regimes of L. In the limit of 
L — y eo, the Kondo temperature is determined by as 
T Koo ~ V T ooU/2 expbr(e d - e F )(U + e d - e F )/(2T 00 U% 
and the cloud size is £,Koo = hv F /(kBTKoo)- For L ^> 
£koc (he., Tk oo 3> A), many resonances of T(e) around 
e F contribute to the Kondo effect, and their average value 
Too results in Tk — Tkoo- In this regime, the Kondo 
cloud is not sensitive to V g . On the other hand, for 
L <C £,Koo (he., Tk oo <C A), the resonance of T(e) located 
at e F , namely To = r(e^), mainly determines the Kondo 
effect, resulting in Tk = Tko ~ \/FoU/2exjp[ir(e d — 
e F ){U + e d - e F )/{2Y Q U)\ and £ K = hv F /(k B T K o). In 
this case, V g directly affects conduction electrons within 
the cloud, and modifies Tk- 

Below, we will discuss the change of Tk between Tko 
and Tkoo as a function of L/^k in the two possible sit- 
uations, case A where one changes L with keeping twD 
(hence Tk oo and £,koo) constant, and case B where one 
changes ty^D with keeping L constant. 

Kondo temperature. — We compute the dependence of 
Tk on L and V g analytically using the poor man scaling, 
and also numerically using the NRG approach [28| . 

In the poor man scaling, the renormalization of J — > 

J + J 2 (J-D + lD°)^ e P( e )/\ e \ i s performed with reduc- 
ing the energy bandwidth of the wire from Dq to D, and 
stopped at the bandwidth where J 2 J dep(e)/\e\ is com- 
parable with J. The final bandwidth provides the scale 
of Tk- From this and Eq. ((3J, we obtain 



T K , eV g cos(k F (2L + a)) 2L 

^ — > — o+ • 2/, \ ^Hf— J> 

±Koa 2t sin (k F aj z,k 



(4) 



eVg -Ci(^) for k F = k Fn , (5) 



2isin (k F a) £,k 
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FIG. 3: (Color Online) Case A under the resonance condition 
of kF — kF, n - In this case, one changes L with keeping twD 
(hence Tkoo and £koo) constant, (a) Kondo temperature Tk 
as a function of L, obtained by the poor man scaling (blue 
solid curve) and NRG (green dashed). The two approaches 
show qualitatively the same overall behavior that Tk drasti- 
cally changes for L < £a-oo, while T K ~ Tkoo for L > £a-oo; 
their discrepancy in L > £,koo is discussed in the text, (b) 
NRG result of the temperature T dependence of conductance 
G for various values of L/^koo- We choose T oa /2t = 0.28, 
e /2t = 0.925, eV g /2t = 0.125, and U/2t = 3.6. 



where Ci(y) = J_^dy'(cosy')/y', k F . n = 2Trn/(2L + a), 
and n is an integer. Equation (U) is obtained by putting 
k -> kp = k(e F ), q ->■ <?(£f), £ > a, v F = i|| | fes , = 
(2at/h) sm(k F a), the small gate voltage of |eVg| <C 
2isin 2 (fc^a), and the linearization of e ~ ep + hvp(k — kp) 
into Eq. ([3]); the simplication of k — > kp and q —y q(ep) 
is valid within the small energy scale of Tk ■ We remark 
that £k depends on L and V g in Eq. (J4j) . 

The dependence of Tk on L in Eq. (j4|) comes from the 
term Ci(2L/£if), which gives information on the cloud, 
and also from the 2kp oscillation by the fact that the 
resonance centers in T(e) shift across tp as L changes. 
One can focus on the former in the situations where the 
latter can be ignored. In case A, where one changes L, 
one can reduce the effect of the 2kp oscillation, by tuning 
the system to a situation that tp is located near the 
bottom of an energy band (near the van Hove singularity) 
where the 2k F term slowly oscillates, or by considering 
the resonance condition of kp = kp^ n where the 2kp term 
provides the maximum value; see Eq. ([5]) . The resonance 
condition can be achieved at each value of L, by tuning 
an additional gate voltage applied to the entire region 
of the wires (not shown in Fig. []} with monitoring the 
conductance through the dot. On the other hand, in case 
B, the term cos(kp(2L + a)) is constant, hence can be 
ignored, since one changes twD with keeping L constant. 

In case A under the resonance situation of kp = 
the poor man scaling in Eq. ([5]) is plotted as a function of 
L/^Koc (rather than L/£k) in Fig.[3Ja). As expected, Tk 
stays at Tkoo for Cifoo ^ L, starts to drastically change 
around L = £koo, an d approaches to Tko for L <C £,k- 
This overall behavior is accompanied by the oscillation 
of \h(Tk /Tkoo) ~ sinc(2L/^ifoo) for L > ^Koo, which is 
obtained from the asymptotic form of Ci(x) ~ (sinx)/x 
for x ^> 1 and by putting £k — £,Koo (valid for L > 



£ifoo) in Eq. (|4]). The sine behavior originates from the 
average effect of the oscillation of p(e) within Tkoo', the 
average becomes more significant with increasing L, as 
the number Tr^/A of resonances within Tkoo increases. 
On the other hand, for L <C £k<x, we find \MTk /Tkoo) °c 
- \n{L/£, K ), using Ci(x) - ln(x) + 0.577 for x < 1. The 
above behavior of Tk(L/£kcc) reveals the Kondo cloud. 

Conductance. — We compute the temperature T de- 
pendence of electron conductance G between the wires 
through the dot, by using the NRG approach 
[28j . Below, we will discuss how to extract £kol 
G(T) in cases A (see Fig.© and B (Fig. 3). 

We continue to discuss case A under the resonance 
condition of kp = kp >n . In Fig. [S^b), we plot G{T) for 



2G 



from 



different Us. It is custom [32| to get an estimate for Tk 
from the temperature at which G(T) equals the half of 
the zero-temperature conductance of the L — > oo case. 
G(T) shows the behavior distinct between L > £koo an d 
L < 6foo- For L > 6foo, G{T) equals to half of G(T = 0) 
of the L — > oo case at almost the same temperature, 
implying that Tk equals Tkoo independent of L in this 
regime. On the other hand, for L < ^oo! G(T) shows 
that Tk changes toward Tko as L decreases. This NRG 
result agrees with the poorman scaling; see Fig. ©a). In 
this way, one directly measures £,Koo . 

There is a discrepancy between the two curves in 
Fig. ©a). In the NRG case, Ixi(Tk/Tkoo) decreases 
only monotonously for L > without showing the 
sine behavior sinc(2L/^R- oc ) of the poor man scaling. 
The discrepancy may come from the known limitation 
that the logarithmic discretization scheme of NRG does 
not perfectly capture the behavior of high-energy states 
(higher than A in our case). On the other hand, for 
L <C £ifoo , where low-energy states mainly contribute to 
the Kondo effect, the NRG shows the same behavior of 
\tl{T k /Tk oo) ~ ~ln(L/^x) as the poorman scaling. 

Next, we discuss case B where one changes twD with 
keeping L constant (hence A = irhvp/L does not alter). 
In Fig. |U we plot the NRG result of G(T) for different 
£ic's. In this case, one can obtain Tk from the high- 
temperature behavior of G{T) in the same way as above, 
by choosing the temperature at which G(T) equals the 
half of G(T = 0) of the L — >• oo case. The result of 
Tk agrees with case A; see the inset of Fig. @|b). We 
below suggest another way to see the cloud, based on 
the low-temperature behavior of G(T). We note that 
as T changes across A, G(T) can show a jump due to 
the resonance structure of p(e), as shown in the case of 
L/6f = 26 in Fig. Ufa). 

We describe the regime of T -c Tk, A, using the fixed- 
point Hamiltonian of the Fermi liquid theory 

Et 1 i e k 
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FIG. 4: (Color Online) Case B. In this case, one changes two 
(hence £k) with keeping L constant, (a) The NRG result of 
G(T) for different values of twD', we show the values of £k 
instead of twD- (b) T c a [defined in Eq. ©] as a function 
of L/£k, obtained from the two different approaches of the 
Fermi liquid theory and the NRG. Inset: The NRG result 
of Tk(L/£koo) exhibits the same behavior as Fig. EJa). We 
choose L = 100a, e /2t = 0.931, eV a /2t = 0.125, and U/2t = 
3.6. 6f is chosen to be located at the center of a resonance of 
r(e), for simplicity. 



where cJ. CT creates an electron with momentum k, spin a, 
and energy e^ , and S p is the phase shift coming from the 
potential scattering (i.e. from the particle-hole symmetry 
breaking). Although p depends on e, we take for simplic- 
ity a constant p in Eq. (|6|) as a crude approximation. The 
second term of Eq. (O describes elastic scattering of elec- 
trons by the Kondo singlet, with scattering phase shift 
5(e) = 7r/2 + (e- ep)/T K + S p . The third term shows 
repulsive interactions between quasiparticlcs that break 
the Kondo singlet, and contributes to inelastic T-matrix 
i in as -7rp(e F )Imt in (e) = [(e - e F ) 2 + n 2 T 2 ]/(2T K ). By 
combining 5(e) and Imi m (e), one can obtain G(T), fol- 
lowing Ref. [fj, as 



G(T) = 



2e 2 1 + cos(25 p ) ir 2 T 2 



2/3 



T- 

2 cos(2<5 p ) 



T 2 



3(a + /3)A 2 (l+cos(2<5 p ))T 2 . 



(7) 

(8) 



where a and (3 are constants depending on V g and kp but 
independent of L. 5 P is obtained by comparing G(T = 0) 
with Eq. (0. When 5 p is nonzero (i.e., when the particle- 
hole symmetry is broken), the potential scattering con- 
tributes to G(T), causing the deviation of G(T = 0) from 
the unitary- limit value of 2e 2 /h. Note that when p is 
independent of e and the particle-hole symmetry is pre- 
served, S p = and (3 = 0, hence, T e g — > Tk- 

T c ff is obtained by comparing G(T) with Eq. ([7]) in ex- 
periments or in the NRG approach, while it is directly 
computed from Eq. © in the Fermi liquid theory. We 
plot T e ff(L/6c) in Fig. gib), which shows that the NRG 
and the Fermi liquid theory are in good agreement; their 
quantitative discrepancy may come from the approxima- 
tion in our Fermi liquid theory mentioned above. 

The dependence of T ff on or L/^Koo is use- 

ful for identifying £k in case B, since A is constant 



so that T c s(L/^k) directly provides the information of 
Tk(L/£k)', see Eq. flS). For L > T g is almost con- 
stant, implying that Tk and £k are independent of L. 
For L < £ K , T cf f (hence T K ) depends on L/£ K . The 
crossover occurs around £r- ~ L, providing a direct evi- 
dence of the cloud. 

Conclusion. — In this Letter we proposed a way to di- 
rectly detect, by electrical means, the Kondo cloud, based 
on the competition between the cloud length and the 
distance between the Kondo impurity and the region of 
conduction electrons to which a gate voltage is applied. 

Our proposal may be within experimental reach. Case 
A, where L varies, may be achieved with keyboard-type 
gate voltages, while one easily tunes two in case B. For 
both the cases, the wire whose Fermi level is located near 
the bottom (van Hove singularity) of an energy band is 
useful to achieve a conclusive evidence of the cloud. Near 
the van Hove singularity, p(e) is sensitive to V g , resulting 
in large difference between T^oo and Tko- Moreover, 
near the van Hove singularity, the Fermi wave length 
1/kp 3> L, £k] hence we expect that the 2k p term in 
Eq. will oscillate slowly as a function of L in the 
regime of L where the transition between Tkoo and Tko 
occurs. If we assume that in this limit 1/kp is also longer 
than the length scale over which V g (x) spatially changes 
at x = L, we expect that then the 2kp oscillations will 
not be washed out. A good candidate for our proposal 
may be a carbon nanotube, where one expects Tk ~ 1 K 
and £k of a micrometer 34|. 
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